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ABSTRACT 

We present a two-component gasdynamic model of an interstellar cloud embedded 
in a hot plasma. It is assumed that the cloud consists of atomic hydrogen gas, inter- 
stellar plasma is quasineutral. Hydrogen atoms and plasma protons interact through 
a charge exchange process. Magnetic helds and radiative processes are ignored in the 
model. The influence of heat conduction within plasma on the interaction between 
a cloud and plasma is studied. We consider the extreme case and assume that hot 
plasma electrons instantly heat the plasma in the interaction region and that plasma 
flow can be described as isothermal. Using the two-component model of the interaction 
of cold neutral cloud and hot plasma, we estimate the lifetime of interstellar clouds. 
We focus on the clouds typical for the cluster of local interstellar clouds embedded 
in the hot Local Bubble and give an estimate of the lifetime of the Local interstellar 
cloud where the Sun currently travels. 

The charge transfer between highly charged plasma ions and neutral atoms gener- 
ates X-ray emission. We assume typical abundance of heavy ions for the Local Bubble 
plasma and estimate the X-ray emissivity due to charge exchange from the interface 
between cold neutral cloud and hot plasma. Our results show that charge exchange 
X-ray emission from the neutral-plasma interfaces can be a non-negligible fraction of 
the observed X-ray emission. 
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1 INTRODUCTION 

During last decades, after the first disc overy of the X-ray 
emission from comets l|Lisse et al. II 19961 ). it was found that 
the X-ray emission produced in charge exchange reactions of 
highly charged ions with neutrals can be observed almost in 
every astrophysical object where the mixing of ions and neu- 
trals takes place. In the solar system X-ray emission has been 
obser ved from many planets - the Earth, Jupiter, Venus and 
Mars. lCravensI (2000) suggested that one of the mechanisms 
of X-rays from these objects is the charge exchange process, 
which highly charged ions of the solar wind undergo in the 
upper atmosphere of the planets. Also it was found that 
the charge exchange of solar wind ions and interstellar and 
geocoronal neutral atoms generates X-ray emission in the 
heliosphere l|Coxlll998r ). lLallement I (|2004f ) has pointed out 
that charge exchange process is the most probable source of 
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X-rays from more distant objects outside of the solar system. 
Since the X-ray emission requires the presence of both neu- 
tral and ionized gases the following phenomena have been 
considered as possible sources of X-rays: (a) galactic wind in- 
teracting with the halo dense cloud; (b) high-velocity cloud 
moving through halo; (c) dense interste llar clouds movin 
in intra-cluster gas. Estimates made by lLallement I (|200 
showed that charge exchange X-ray emissivity from spatially 
small regions of the neutral-plasma interfaces may have the 
same order of magnitude as the X-ray emissivity from the 
hot gas itself. 

In order to use X-rays produced by charge transfer 
for determining the physical properties of plasma-neutral 
interfaces, these regions need to be studied theoretically. 
Moreover, the charge exchange changes the structure of 
the plasma-neutral interfaces. The region where dynamical 
effects of the charge exchange are intensively explored is 
the heliospheric interface, i.e. the interface where the solar 
wind interacts with the local interstellar medium. Due to 
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the charge exchange with both interstellar and solar wind 
protons the interstellar neutrals provid e additional pressure 
movi ng the structure closer to the Sun ()Baranov fc Malamal 
Il993h . The resonant charge exchange is the driver for many 
physical phenomena i n the outer heliosphere and at the he- 
liospheric boundaries (jlzmodenov et aill2009l ). 

Often in exploring the physics of the neutral-plasma in- 
terfaces it is necessary to employ a two-component approach 
since the neutral and charged components are not in equi- 
librium. Also local equilibrium may or may not occur inside 
each of the components. If the component is in local equi- 
librium then it can be described as a fluid, otherwise a ki- 
netic approach is required. For example, a kinetic approach 
is needed to describe the interstellar neutral component in 
the heliosphere since the mean free path of the neutrals is 
comparable with the siz e of the heliospheric interface region 
l|lzmodenov et alj|200d ). 

The goal of this paper is to make an initial step in theo- 
retical expl oration of the astr ophysical phenomena that were 
selected bv lLallement I (|2004f ) as potential sources of charge- 
exchange induced X-ray emission. Doing this we develop a 
two-component model of cold neutral cloud embedded into 
a hot plasma. In this initial study we explore the role of 
proton-neutral charge exchange on the formation and time 
evolution of neutral-plasma interfaces. 

The models of cold cloud - hot plasma interaction have 
been developed for decades and the main goal of these mod- 
els wa s to estimate the lif e time o f clouds ( Graha m fc Langen 
ill973t). ICowie fc McKee! (|l977l l. iNagashima et all (|2006l )7 
I Vieser fc Henslerl 1 2007 )). The major difference of our 



model from the existing models is the two-component ap- 
proach that we employ as compared with the one-fluid ap- 
proach used before. 

The pioneering theoretical models have consid- 
ered a cold neutral c loud embe dded in a warm gas 
Zel'dovich fc Pickel'ne 3 l| 19691 ). IPenston fc Brownl 



1973)) and assumed pressure equilibrium and thermal 



balance of the cold and warm gases. The model takes into 
account the thermal conduction and radiative processes of 
heating and cooling of interstellar medium. The condition 
when steady solution exists was found. This model has bee n 
applied for spherical clouds in iGraham fc Langerl l|l973h . 
They considered the clouds with the temperature T ~ 100 
K and the number density n ~ 100 cm -3 . Surrounding 
warm medium has the temperature T ~ 10 4 K and the 
number density n ~ 0.1 cm -3 . For these conditions they 
gave the evaporation rate as a function of a cloud radius 
and medium pressure which allows to estimate the timescale 
on which the clouds evaporate. Small clouds with radius 
0.02 < Rdoud < 0.03 pc at the pressure p/k = 10 3 Kcm~ 3 
evaporate in ~ 10 7 years. 

Later it was found that cold neutral clouds are sur- 
rounded by very hot ionized gas. Because of the fact that 
electrons have high thermal v elocities the heat condu ction 
should be taken into account. I Cowie fc" McKed l| 19771 ) have 
considered classical and saturated heat conduction. In ab- 
sence of magnetic fields and radiative processes the evap- 
oration rate of sph erical clouds was derive d for tw o cases . 
Subsequent work s (iMcKee fc Cowiel (Il977f ) Isiavml (|l989h . 
iDalton fc Balbusl (|l993h . l Vieser fc HensleU <|2007t >> have ex- 



tended these models for two and three dimensional geometry 



and applied magnetic fields, radiative heating and cooling 
processes and different types of heat conduction. 

The charge exchange process (H + H + < — ► H + + H) 
was not considered in either of these models. That was be- 
cause both components have been treated as being in local 
equilibrium or in a state close to the local equilibrium. In this 
case the charge exchange does not have any pronounced ef- 
fect. The situation is different when the components are not 
in equilibrium as in the case of the heliospheric interface. 

When neutral hydrogen gas and hydrogen plasma en- 
counter one another the main process that starts to act is 
the charge exchange process. This is because the charge- 
exchange cross section is larger as compared with the corre- 
sponding (i.e. momentum transfer) cross section for elastic 
collisions between the components. The charge exchange is 
the most effective for momentum exchange between the com- 
ponents. Since the momentum exchange term in the momen- 
tum equation for each of the components is proportional to 
the velocity difference between the components (multiplied 
by the cross section and the number densities of both com- 
ponents) this term acts as a drag and in fact should reduce 
the velocity of the neutral component. This may results in 
the increasing of the cloud lifetime as compared with the 
models where the charge exchange is neglected. 

Another important feature of the model proposed here 
is the assumption that although the total pressure in the 
cold cloud is assumed to be in balance with pressure in hot 
plasma (P c ioud = Ppiasma) there is no pressure balance for 
each of the components. If we ignore the charge exchange 
there is nothing that would prevent the neutral cloud from 
expanding to the ambient plasma and plasma from filling 
the cloud cavity. As we will see below the charge exchange 
is a very efficient process preventing the components from 
fast penetration to each other. 

In Sect. 2 we introduce the two-component model of 
the interaction of cold neutral cloud and surrounding hot 
plasma. The problem is formulated in two cases - "adiabatic" 
and isothermal. A gasdynamic structure of the plasma- 
neutral interface is discussed in Sect. 3. Estimates of charge 
exchange X-ray emission from the interaction interface are 
given in Sect. 4. The lifetime of the interstellar clouds in- 
side the Local Bubble obtained in the frame of the two- 
component model is discussed in Sect. 5. Conclusions are 
drawn in Sect. 6. 



2 MODEL 

2.1 Adiabatic model 

We consider a problem of the interaction of a cold neutral 
cloud and hot fully-ionized plasma. The cloud is assumed to 
be spherically-symmetric and consists of atomic hydrogen 
gas. Plasma is quasi-neutral and consists of protons, elec- 
trons and highly charged ions (n e ~ n p + ni). We assume 
that rii/np <C 1. The medium of hydrogen atoms and ions is 
described in the two-fluid approach. H atoms and plasma 
ions interact solely through the charge exchange. In this 
subsection we present an "adiabatic" model of the cloud - 
plasma interaction. It implies that heat fluxes in the plasma 
component are ignored, except for the energy transfer be- 
tween two components due to the charge exchange. 



Two- component model of the interaction of an interstellar cloud with surrounding hot plasma 3 



Fluid approach for H atoms is reasonable since (as it is 
seen below) the frequency of the collisions between H atoms 
exceeds the frequency of the charge exchange between H 
atoms and plasma protons almost in the entire interaction 
region. This takes place due to relatively large number den- 
sity of neutrals. Coupling of the plasma and neutral compo- 
nents by the charge exchange results in the momentum and 
energy exchange between plasma and neutral gas which are 
not in local thermodynamic equilibrium between each other. 
In the model we assume that effective collisions between H 
atoms lead to local thermodynamic equilibrium in the neu- 
tral component and local thermodynamic equilibrium in the 
plasma is supported by "collective" processes. 

The governing equations for plasma and hydrogen gas 
flows are the following: 

dm 1 dr 2 (nu)i _ ^ 
dt r 2 dr 
d{nu) l 1 d(r 2 (nu 2 +p)i) 2 Pl 

— al 1" ~2 a = 1~ 9 2 < ' 

at r z or r 

d(n(e + 1 9 (r 2 (nu(e + ^) + pu)?j ^ 

dt r 2 dr ' 

where i = p for the plasma and i — H for the hydro- 
gen gas. Here n p and tih are the number densities of pro- 
tons and H atoms, respectively; e p = p P [(7 — and 
en = PhKj — 1)ph]~ are specific internal energies of plasma 
and hydrogen gas. It is assumed that in the plasma the elec- 
tron and proton temperatures are equal. The plasma and hy- 
drogen gas temperatures can be determined from the equa- 
tions of state p p = 2ripkTp and pu = nnkTu, respectively. 

The terms , q3 t are the sources of momentum and 
energy due to the charge exchange of H atoms and plasma 
protons. For the so urce terms we use the expressions by 
iMcNutt et all (|l998h : 

q2 p = -q2 H = n p v H {uH - u p ) 

q Sp = -qs H = npvn + Jj^(T H - T p ) j 

Here uh denotes the charge exchange frequency: uh = 
riHU<jfj'{U), where a^ x p (U) is the charge exchange cross 
section and U is the relative velocity of H atom and pro- 
ton. Fo r charge exchange cross section we use the expres- 
sion bv lMaher fc Tinslevl (|l977l ): a^ p = (1.64 ■ 1(T 7 - 6.95 ■ 
10~ 9 lnU) 2 . The exp ressions for U* and U can be found in 
IMcNutt et all l|l998T l. 

The initial conditions are the following: 

< r < R c : n H = n H1 ,T H = Ti,uh = 0,n p = 

r > R c : rip — n P 2,Tp = T2, u p = 0, ns = 

At the initial moment of time a spherical cloud with the 
radius R c is filled by the cold hydrogen gas (Th ~ 100 — 10 4 
K) and surrounded by the hot plasma (Tp ~ 10 6 K ). It is 
assumed that at t = the cold neutral cloud and the hot 
plasma are in pressure equilibrium: nniTi = 2n P 272. Since 
Ti < T 2 then n H i < n p2 . 

We solve the problem in dimensionless form. For the 
characteristic length the mean free path of H atom in hot 
plasma is adopted: L = [no<r^ ; p (?io)]~ 1 . For the character- 
istic velocity it is convenient to take the thermal velocity 



of hot protons: no = ^2kT 2 /m p . The number densities are 
normalized to the proton number density in hot plasma: 
no = n P 2 . Then it is easy to show that the solution depends 
only on two dimensionless parameters: n = nHi/n P 2 and 

Rc — Rc/ L- 

For the simplicity of the numerical solution of the prob- 
lem we introduce small portions of protons and H atoms in 
the cloud and hot gas respectively. n pl is the number den- 
sity of plasma portion in the neutral cloud and riH2 is the 
number density of H atoms in surrounding plasma. These 
values are close to zero: n p i <C rim, riH2 <^ n P 2 and do not 
affect the solution. 

The problem is solved numerically by the Godunov 
shock capturing method with improved order of scheme ac- 
curacy. The space grid has been reconstructed on each time 
step according to the propagation of a shock wave, a contact 
discontinuity and a rarefaction wave-front in either neutral 
or plasma component. To improve the scheme resolution we 
used the flux limiter function "min-mod" . 

2.2 Isothermal model 

For high temperature plasma the electron heat conduc- 
tion plays an important role since hot electrons with large 
thermal velocities penetrate into the cooler regions of the 
plasma and heat it by collisions. The thermal conductiv- 
ity in a fully ionized hydrogen plasma is given by l|Spitzerl 
Il957l ) :\c = Xc(T p ) = aTp /2 ergscm~ 1 s~ 1 K~ 1 , where 
a ~ 10~ ergsan' -^s^K 112 . One can show that since the 
temperature of the Local Bubble plasma is about 10 6 K 
then the heat flux q — —XedT p /dr is the order of magnitude 
above the other terms in the energy equation for the plasma. 
Therefore, we can consider the extreme case, when the elec- 
tron thermal conduction leads to instant plasma heating in 
the cloud-plasma interface, and the plasma flow can be de- 
scribed as isothermal. 

The system of equations and initial conditions govern- 
ing the cold cloud interaction with hot isothermal plasma 
has the same form as for "adiabatic" model, except the en- 
ergy equation for the plasma component which is replaced 
by the condition: T p (x,t) = const. 

3 GASDYNAMIC STRUCTURE OF THE 
NEUTRAL CLOUD-PLASMA INTERFACE 

3.1 Adiabatic model 

In our calculations we varied the dimensionless parameters 
n and R c . Here to discuss a qualitative gasdynamic struc- 
ture of a cloud-plasma interface we present a solution for 
the values n — 285, 1143 and R c = 20. It corresponds to 
all situations when the number densities in cold neutral 
cloud and surrounding hot plasma satisfy the conditions 
n_sr=285n p and xih = 1143n p , respectively. The parameter 
n and pressure balance determine the ratio of temperatures 
in the cloud and the plasma Sf 1 . For example, if we as- 

sume n p = 0.0009 cm' 3 ,T p = 10 6 K for the Local Bub- 
ble plasma then h = 285 corresponds to the clouds with 
n H = 0.25 cm,- :i ,T H = 7000 K and n = 1143 - to the 
clouds with rtH = 1 cm -3 ,Th = 1750 K. Radius of the 
cloud is 2.9 pc. 
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Figure 1. Distributions of dimensionless number density, velocity, pressure and temperature of neutral gas (top plots: A-D) and plasma 
(bottom plots: E-H) as functions of a radial distance at t = 1,5. Solid curves correspond to the solution for h = 285, dashed curves - 
n = 1143. R c = 20. (S) and (CD) denote a shock wave and a contact discontinuity, respectively. 




After the beginning of the interaction the neutral gas 
flows into the hot plasma and the plasma penetrates into the 
neutral cloud. The coupling of the plasma and neutral gas 
in the charge exchange pr ocess substantially influe nces the 
flows of both components l|Provornikova et al.ll2009l ). Figure 
[T]shows number densities, velocities, pressures and tempera- 
tures of the hydrogen gas (plots A-D) and the plasma (plots 
E-H) as functions of a radial distance for two moments of 
time. After the beginning of the interaction a contact discon- 
tinuity (CD) and a shock wave (S) are formed in the neutral 
gas and propagate into the hot plasma. It can be seen from 
plots A-D for t — 1 in Fig. [1] The neutral gas is compressed 
and heated behind the shock. Distributions for t = 5 show 
that the charge exchange of H atoms and plasma protons 
leads to that the shock strength decreases, and gasdynamic 
parameters of the neutral gas change continuously. 



In the plasma a shock wave (S) and a contact discon- 
tinuity (CD) propagate into the neutral cloud at the initial 
stages of the interaction (Fig. [1] E). The discontinuities are 
decelerated very fast due to the interaction with neutral gas 
by the charge exchange. The plot E in Fig. [1] shows that at 
the edge of the cloud the plasma density has a maximum 
several times larger than the number density in surrounding 
hot plasma. Plots B, F in Fig. [T]show that due to the ex- 
change of momentum between the components in the charge 
exchange process the velocities of neutral gas and plasma de- 
crease with time. It is important to note that at the edge of 
the cloud the narrow region is formed where the tempera- 
ture of both components gradually changes from its value in 
the cold cloud to the value in the hot plasma. Distributions 
of hydrogen gas parameters obtained for different values of 



Two- component model of the interaction of an interstellar cloud with surrounding hot plasma 5 




Figure 3. Distributions of gasdynamic parameters of neutral gas (top plots: A-D) and plasma (bottom plots: E-H) as functions of a 
radial distance for t = 1,5. Solid curves correspond to the solution in the isothermal model, dashed curves denote the solution in the 
"adiabatic" model. Here h = 285 and R c = 20. (S) and (CD) denote a shock wave and contact discontinuity, respectively. (Si) shows a 
shock wave formed in plasma. 




Figure 4. Space-time distributions of number density, velocity and temperature of neutral gas (top plots: A-C) and number density, 
velocity of plasma (bottom plots: D-E). Plot F shows a charge exchange rate of H atoms and plasma protons. The solution corresponds 
to n = 285 and R c = 10. 



the parameter h show that the solution for the neutral gas 
beyond the cloud does not depend on the parameter h. 

Figure[2]presents space-time distributions of the neutral 
gas and plasma parameters. The results show that the inter- 
action of the neutral gas and hot plasma by charge exchange 
leads to the formation of the interface region at the edge of 
the cloud. From the simulation up to t — 1000, which cor- 



responds to the 1.1 Myr, it is seen that the interface exists 
during rather long time periods. The neutral cloud-plasma 
interface is characterized by the following features: 

1) plasma number density has a maximum at the cloud 
boundary (Fig[2] E). In the interaction region the plasma 
moves towards the cloud boundary and decelerates due to 
momentum exchange with H atoms ( Fig. [2] F). It results 
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in increasing of plasma density with time. For instance, over 
the interaction period of 1.1 Myr plasma number density 
increases by 15 times in the transition region. 

Since the number density of the hydrogen gas in this 
region is the same order as the density in the cold cloud 
(Fig. El A) the frequency of charge exchange between plasma 
ions and H atoms has a maximum here; 

2) velocity of neutral gas in plot B of Fig. [2] shows that 
the cloud does not expand during the interaction with sur- 
rounding plasma; 

3) temperature of neutral gas shows that the cloud does 
not suffer from heating (Fig. [21 D) and even at large time 
intervals the gas inside the cloud has an initial tempera- 
ture Th ■ The charge exchange process results in the same 
temperature of neutral and plasma components throughout 
the interface. The temperature increases from cold tempera- 
ture of the cloud Th to the high temperature of hot plasma 
T p (Fig. [J] D,H). Numerical solution shows the existence of 
waves which periodically appear at the edge of the cloud, 
propagate into the cloud and than reflect back from the 
cloud center. 

Thus in the "adiabatic" model the cloud remains cold 
and its boundary is not displaced over a long period of time. 
It will be seen below that taking into account a heat conduc- 
tion in the plasma strongly affects this gasdynamic structure 
of the cloud-plasma interface. But from results of the "adi- 
abatic" model one can conclude that the charge exchange 
process may be considered as the most important mecha- 
nism providing the existence of the interstellar clouds in the 
Local Bubble hot plasma. 

3.2 Isothermal model 

The isothermal model implies that the plasma flow is 
isothermal everywhere including the interaction region: 
T v (r,i) = T%. Figure [3] shows the distributions of gasdy- 
namic parameters for neutral gas and plasma as functions 
of a radial distance for two moments of time. In fig. [3] we 
also present the results of the "adiabatic" model for the sake 
of comparison. 

When the cloud gas starts interacting with the ambient 
plasma the structure of the neutral and plasma flows is sim- 
ilar to the "adiabatic" case. A contact discontinuity (CD) 
and a shock (S) are formed in the neutral gas and propa- 
gate into the hot plasma (Fig. [3l A) . Since the plasma flow is 
isothermal only a shock wave exists and tends to move into 
the cloud (Fig. [3] G). Plot A in Fig.[3]shows that the number 
density of H atoms increases at the edge of the cloud and 
velocity becomes negative(Fig.[3] B). Since H atoms interact 
with the hot protons by charge exchange the temperature of 
H atoms increases in the interface. The narrow region where 
H atoms are heated and have a maximum density propa- 
gates inside the cloud. On the other hand due to the heat 
flux into the neutral component from hot plasma the pres- 
sure of neutral gas increases in the interaction region and 
gas flows out into the hot plasma. 

Momentum exchange of plasma protons and H atoms 
results in that the plasma velocity becomes positive almost 
in entire interaction region (comparing to the "adiabatic" 
case where in the interface plasma flows towards the cloud 
with decreasing velocity) (see Fig. [3] F). Thus the neutral 
gas expanding to the ambient plasma captures the plasma 



flow. The region at the edge of the cloud where the plasma 
density and velocity has a maximum (Fig. [3] E, F) propa- 
gates to the hot medium. Steepening of that region creates 
a shock wave in the plasma (Si) (Fig. [3] E). 

Comparison of the solutions for the isothermal and "adi- 
abatic" models in fig[3] also shows that in the isothermal 
model disturbances in neutral gas propagate faster than in 
the "adiabatic" model. That occurs due to the heat transfer 
from hot isothermal plasma into the neutral gas. 

Figure [4] shows space-time distributions of the neutral 
gas and plasma parameters. In calculations n — 285, R c = 
10. One can see in figure [4j plot A that in the isothermal 
model the size of the cloud changes with time. The cloud 
expands due to the heat flux from hot isothermal plasma, 
then cloud radius decreases and the cloud vanishes. 

The main results of the modeling of the interaction of 
cold neutral cloud and hot isothermal plasma are the follow- 
ing: 

1) neutral cloud is being heated with time due to cre- 
ation of new hot H atoms from hot protons in the charge ex- 
change process. Space-time distributions of neutral gas tem- 
perature shows that the temperature of the cloud is equal 
to the temperature of surrounding plasma at t ~ 1300 (Fig. 
®C); 

2) inside the cloud the density of neutral gas decreases 
with time but plasma density increases (Fig. 2] A, D). Over 
the period t w 1300 number density uh changes from the 
initial value uh = 285 to hu = 1 in the cloud, number 
density of plasma reaches h v = 0.1; 

3) The velocities of both neutral and plasma compo- 
nents decrease with time due to momentum exchange be- 
tween the components. 



4 ESTIMATE OF CHARGE EXCHANGE 

X-RAY EMISSION FROM CLOUD-PLASMA 
INTERFACE 

The distributions of the neutral gas and plasma parameters 
obtained in the frame of the models allow us to calculate 
the rate of charge exchange between H atoms and plasma 
protons f3 = nHn p Ua^ x p for the plasma-cloud interaction re- 
gion. Figure [S] shows a space-time distribution of the charge 
exchange rate obtained in the "adiabatic" model. The charge 
exchange rate has a maximum inside the cloud-plasma in- 
teraction interface. This maximum exists due to increased 
number densities of both plasma and H atoms in this re- 
gion. This increase of charge-exchange rate is an indication 
that the charge exchange X-ray emission from the neutral- 
plasma interface is not negligible as it was pointed out by 
lLallement I (|2004h . 

Let us estimate an importance of the charge exchange 
X-ray emission from the interaction interfaces as compared 
with the hot gas X-ray emission. The charge exchange X-ray 
emissivity is expressed by the integral: P ex — J^ 2 ftidr = 

J^ 2 riinHUa^Jdr, where r\,V2 are boundaries of the tran- 
sition region, /3j is the rate of charge exchange between H 
atoms and plasma heavy ions, m is the number density of 
heavy ions, cr^. 1 is the charge exchange cross section of H 
atoms and heavy ions. To get a simple estimate we assume 
that X-ray photons are generated in the charge exchange 
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Figure 5. a) Space-time distribution of charge exchange rate 
H — H + in "adiabatic" model, b) Charge exchange rate as a 
function of a radial distance at t = 1000. Charge exchange 
rate has a maximum in the transition region for both "adia- 
batic" (solid curve) and isothermal (dashed curve) models. Here 
n = 285, R c = 20. 

process between H atoms and oxygen ions +r and an abun- 
dance of ions in the hot plasma is : n +7 /n p ~ 6 ■ 1CP 5 . 
Here the typical abundance of oxygen ions +7 for solar 
wind is given. The ratio of charge exchange cross sections is: 
a?° +7 ~ 1.7 l|Koutroumpa et alJl200d ). 
Therefore 

_ n 0+ r o ex 

^ ex ~ Z Hp lcx ' 

where I ex = J^ 2 n p nHUcr ex p dr is the integral of charge ex- 
change rate of H atoms and plasma protons which can be 
calculated by using our model results. Figure[6]shows the de- 
pendence of I ex on time. In the "adiabatic" case for the Local 
interstellar cloud with parameters tin = 0.25 cm -3 ,Th = 
7000 K corresponding to dimensionless parameter h = 285 
one can estimate: 

P ex « 10 3 keVcm^s' 1 . 

The volume emissivity of diffuse X-ray emission from hot 
gas can be estimated as: Ph = 5.8 ■ W~ 14 n 2 keVcm~ 3 s _1 , 
where n e is the n umber density of electrons in hot plasma 
ijLallement 1120041 1. In our calculations n e — 0.0009 cm -3 . 
Then Ph ~ 4 ■ 10 -20 keVcrn~ 3 s -1 . Assuming the character- 
istic size of the Local Bubble is L lb ~ 100 pc the emis- 
sivity of diffuse X-ray from the hot plasma is PhLlb ~ 
WkeVcm" 2 s" 1 . The emissivity ratio is: P ex / 'PhLlb ~ 10 2 . 
The charge exchange X-ray emissivity from the cloud- 
plasma interface is two orders above the thermal emissivity 
from hot gas. Since "adiabatic" model does not take into 
account thermal conduction processes the calculated mag- 
nitude of charge exchange X-ray emissivity is overestimated. 

Results of the isothermal model also show that the 
charge exchange rate has a maximum in the neutral cloud - 
plasma interface (Fig. U F) , but the value of the maximum 
is much less than in the "adiabatic" model (Fig. [5] b) . Using 
the average value of I ex (see fig. [6] dashed curve) for the Lo- 
cal interstellar cloud the charge exchange X-ray emissivity is 
estimated as: P ex « 10 keVcm~ 2 s _1 . Using the estimate for 
the diffuse hot plasma X-ray emissivity from the Local Bub- 
ble one can see the ratio P ex / PhLlb ~ 1. Thus even in the 
isothermal case (which provides us low limit of X-rays) the 
emissivity from the neutral-plasma interface due to charge 
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Figure 6. Time variation of an integral of H — H + charge 
exchange rate. The solid curve corresponds to the "adiabatic" 
model, dashed curve - isothermal model, n = 285. 

exchange has the same order of magnitude as compared with 
the hot gas X-ray emissivity. 



5 LIFETIME OF THE INTERSTELLAR 
CLOUDS 

Let us define a lifetime of a cloud as a timescale over which 
the temperature in the cloud becomes equal to the tem- 
perature in surrounding hot medium. Our two-component 
model allows to give an estimate of the lifetimes of inter- 
stellar clouds surrounded by hot Local Bubble plasma. This 
can be done for two extreme cases - the "adiabatic" and 
isothermal. 

In section 2 we showed that in the frame of the "adia- 
batic" model the cloud boundary is nearly steady and the 
cloud may exist more than several Myr. For example, simu- 
lation of the Local interstellar cloud (LIC) with parameters 
n H = 0.25 cm~ 3 ,T H = 7000 K embedded in the Local 
Bubble with n p = 0.0009 cm~ 3 ,T p = 10 6 K shows that 
the LIC lifetime may reach more than 5.6 Myr. This esti- 
mate changes essentially when an electron heat conduction 
is taken into account. 

We performed model calculations for different radii of 
the cloud. Figure [7] shows the dependence of the cloud life- 
time on the cloud radius. The results are presented for 
two values of the dimensionless parameter h = 285, 1143. 
Solid curves denote the results of two-component isother- 
mal model. The lifetime of the clouds with the number den- 
sity ~ 0.25 cm -3 and temperature ~ 7000 K is estimated 
as 0.2-1.5 Myr. The colder clouds with the number density 
~ 1 cm -3 and temperature ~ 1700 K may exist twice as 
long: the lifetime of the clouds with the radius ~ 1.1-1.5 
pc are ~ 1.6-2.7 Myr. It is known that the Local interstel- 
lar cloud has irregular shape with the characteristic scale 
about 1-10 pc. If we adopt the size of the LIC as 3 pc then 
its lifet ime is estimated a s ~ 1 .5 Myr. 

In[c owie fc McKeei (1 19771 ) the analytical formulae for 
the lifetime of the spherical neutral clouds embedded in hot 
plasma were derived. They considered two types of heat 
conduction in fully ionized hydrogen plasma: 1) classical 
thermal conduction q ~ T 7//2 ; 2) saturated heat conduction 
q ~ T 3/2 , when a mean free path of electrons in plasma is 
comparable to the temperature scale height. Magnetic fields 
and radiation processes were ignored in the model. Dashed 
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Figure 7. The lifetime of interstellar clouds as a function of 
cloud radius. Black curves correspond to the solution for n = 285, 
grey curves - n = 1143. Solid curves denote t he results of two- 
compo nent isothermal model, dashed curves - ICowie fe McKeel 
( 1977) in the case of clas sical thermal conduction, dash-dot curves 
- ICowie fc McKeel (l977h in the case of saturated heat conduction. 

curves in fig. [7] show the lifetime calculated in the case 1). 
Dash-dotted curves denote the results obtained in the case 
2). For the clouds with the radii greater than 1.1 pc the 
lifetime obtained in t he two-component mod el exceeds the 
lifetime calculated in ICowie fc McKeel <| 1977ft in both cases 
1) and 2). For denser clouds (n = 1143) with radii 0.1 — 1.1 
pc the lifetime in the two-component model is comparable 
with the lifetime in the case 2) with the saturated heat con- 
duction. Note, that our two-component model takes into ac- 
count the charge exchange process and the plasma heating 
by the heat conduction. The latter process is considered in 
extreme case of the isothermal plasma flow. For more re- 
alistic situation the effect of heat conduction is expected 
to be somewhat smaller and, therefore, the cloud lifetime 
will be larger. Nevert heless, the comparison o f lifetime esti- 
mates with the models lCowie fc McKeel l|1977l ) demonstrates 
that the lifetime i n our model is larger or comparable with 
ICowie fc McKeel lll977ft . That means that the charge ex- 
change process plays a significant role in the evolution of 
the interstellar clouds and should be taken into account. 



6 CONCLUSIONS 

We have presented a two-component gasdynamic model of 
the interaction of a cold neutral cloud and surrounding hot 
plasma. It was assumed that the neutral and plasma compo- 
nents interact by the charge exchange process. The problem 
has been considered in two extreme cases - "adiabatic" and 
isothermal. In the "adiabatic" model the heat conduction in 
plasma is ignored. In the isothermal model it is assumed that 
the plasma is isothermal due to effective electron heat con- 
duction. We obtained gasdynamic solutions which describe 



the structure of the neutral-plasma interfaces. The results of 
the "adiabatic" model can be briefly summarized as follows: 

• The interaction of cloud neutral gas and plasma 
through charge exchange process results in the formation 
of a plasma-neutral interface, i.e. the transition region at 
the edge of the cloud; 

• In the transition region the number densities of both 
neutral and plasma components increase significantly and 
become non-negligible; 

• Inside the transition region the temperature of the 
medium (for both plasma and neutral) has all of the in- 
termediate values; 

• Over the entire period of the cloud-plasma interaction 
the cloud remains cold and its boundary is nearly steady; 

• The emissivity of X-ray generated by the charge ex- 
change process of plasma heavy ions and H atoms in the 
neutral-plasma interface is two orders of magnitude larger 
than the emissivity of diffuse X-ray from the hot medium. 

The basic results obtained in the isothermal model are 
the following: 

• The interaction of the neutral cloud with hot isothermal 
plasma results in heating and disappearing of the cloud; 

• The lifetime of the interstellar clouds with the number 
density ~ 0.25 cm -3 and temperature ~ 7000 K surrounded 
by the Local Bubble plasma with n p — 0.0009 cm~ i and 
T p = 10 6 K is estimated as 0.2 - 1.5 Myr. The lifetime of 
the Local interstellar cloud is estimated as 1.5 Myr; 

• We conclude also that the charge-exchange X-ray emis- 
sivity from the neutral-plasma interfaces may have the same 
order of magnitude as diffuse X-ray emissivity from the ex- 
tended hot plasma regions. It means that in the analysis of 
experimental data on X-ray emission in different astrophys- 
ical cases where the neutral or partly ionized gas interacts 
with plasma it is important to consider the contribution of 
X-rays induced by the charge exchange. 
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